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Abstract 

A generalised ladder operator is used to construct the conserved 
operators for any one-dimensional lattice model derived from the Yang- 
Baxter equation. As an example, the low order conserved operators 
for the XYh model are calculated explicitly. 

1 Introduction 

The method for constructing integrable one-dimensional lattice models from 
solutions of the Yang-Baxter equation is well known (eg. see [1]). In principle, 
the conserved operators can be obtained by series expansion of the family of 
commuting transfer matrices. A more practical approach is to use the ladder 
operator which permits a recursive method through repeated commutators 
to obtain the conserved operators. 

For models where the solution of the Yang-Baxter equation has the dif- 
ference property, it has been established [2, 3, 4] that the ladder operator is 
a lattice analogue of the boost operator for Lorentz invariant systems. Re- 
cently it has been shown that for the Hubbard model, which is not Lorentz 
invariant in the continuum limit as a consequence of spin-charge separation, 
and is reflected by the fact that the solution of the Yang-Baxter equation 
does not have the difference property, the ladder operator still exists [5]. 

The present work extends [5] to develop a general theory for the construc- 
tion of the ladder operator for any integrable system obtained through the 
Yang-Baxter equation. The theory will be applied to analyse the conserva- 
tion laws for the XY model in a transverse magnetic field. 

2 Integrable Lattice Models using the Quan- 
tum Inverse Scattering Method 

We begin with a vector-dependent solution of the Yang-Baxter equation 
Ri 2 {u, v)R 13 (u, w)R 23 (v, w) = R 23 {v, w)R 13 (u, w)R 12 (u, v) 
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where u, v and w are m-component vectors. Throughout, we assume the 
regularity property R(u, u) = P. Define a set of m local Hamiltonians 



dm 



i = 1, .., m 



with the corresponding global Hamiltonians acting on a one-dimensional lat- 
tice of length L given by 

1=0 

Throughout, periodic boundary conditions are assumed on all summations 
which are evaluated over the length of the lattice. Note it is implicit that all 
the operators h{i} are in fact functions of v. 
The transfer matrix is constructed through 

T(u,v) = tr a (^Ra(L-i) (u, vj...R a ± (u, v)R a Q(u, vj^j 

where a refers to the auxiliary space, which by the standard argument gives 
rise to a commutative family in the first variable; i.e. 

[T(u,v), T{w, v)} = 0, Wu,w. (1) 

It can also be easily verified that 

[H{i},T(u,v)} = 0, Wu. (2) 

It is convenient, however, to define the conserved operators as 

Qni+...+n m 

where they appear in the series expansion 

lnT{u,v)=}^ ■ ■ t{n\. (3) 

V n l \...n m \ 

Thus it follows from (1) that 

[t{n}, t{k}} = 0, Vn,£ 

and moreover from (2) 

[H{i},t{n}] = 0, Vi,n. 

Note that n is an m-component vector with non-negative integer entries. 
Introducing the notation {e i } r [L 1 for the basis of the m-dimensional vector 
space, we can write 

m 

n = J2 n i?i- 
1=1 
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3 Recursion Formula for Calculating the Con- 
served Operators 

For each of the index labels % we define a ladder operator 

1=0 

with the coefficients I taken from the set of integers modulo L. For any 
function </> admitting a Taylor's series expansion we have 

[B{z}, <P(T)} = T.H{i}4\T) 

where T = T(u, u) and 4>' denotes the derivative of (f>. Choosing <fi to be the 
logarithm now gives 

[B{ l },\nT] = H{t}. 

It can be shown that 

[B{i}, T(u,v)} = W 

As a result we obtain the following recursion formula from (4) and the ex- 
pansion (3) 

t{n + e^=[B{i},t{n}] + ^ (5) 
The first few terms in (3) can be identified immediately 

t{6} = \nT, t{$} = H{i}. (6) 

In principle, through repeated use of (4) expressions for all the operators 
t{n} may be obtained. 

Applying the recursion (5), the second order conserved currents can be 
obtained by the following formula: 

*{3 + 5) = \H Mo], hi-!{i}] + \ £[M0, fy-iO'}] 
l dH{j} l dH{i} 

2 dvi 2 dvj ' 1 ' 



3 



4 The XYh Model 



The XY model in a transverse magnetic field has the following Hamiltonian: 

N 

H = ^(Jxofof+i + JyO-i^i+i + ha t) J x, Jy, h const. 

1=1 

This model is known to be integrable [6]. Barouch and Fuchssteiner [7], 
Araki [8], and Grabowski and Mathieu [9] have explicitly calculated the low 
order conserved operators. These results have been reproduced using the 
generalised ladder operator method. 



4.1 R Matrix of the XYh Model 

Bazhanov and Stroganov [10] constructed an elliptic parametrization for the 
Boltzmann vertex weights of the XYh model. In this parametrization, the 
weights are meromorphic functions of 3 complex variables, u — (ui,u 2 ) and 
v = (i>i, 1*2), where only the first vector entry contains the difference property. 
The R matrix is 



R(u, v) 



(R\\ 








^22 \ 





pl2 
K 12 


E>12 
^21 








p21 
K 12 


p21 
^21 





\R\l 








^22/ 



with 

c>22 
-^22 

p21 
^21 



pl2 
^21 



pll 
rl 22 



p(l - e(ui - Vi)e{u 2 )e(v 2 )) 
p(e(u 2 ) - e(ui - vi)e(v 2 )) 

p\j e(u 2 )s{u 2 ) y/e(v2)s(v 2 )(l - e(«i - Vi)) 



R{{ ■■ 

pl2 
K 12 



p{e(ui - uj) - e(u 2 )e(v 2 )) 
p(e(v 2 ) - e(ui - vi)e(u 2 )) 



p21 _ 
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= Rfj = -ikp^J e(u 2 )s(u 2 )^J e(v 2 )s(v 2 )(l + e(«i - Ui))s 



Ui - Vi 



where p is an arbitrary constant and s and e are the respective elliptic func- 
tions sn and (cn + i sn). By imposing the condition R(u, u) = P, we obtain 

the value of p = 1 _ e 2 ( - t)2 - ) ■ 



4.2 Local Hamiltonians 

The local Hamiltonians h t {l} and hi{2} are given as follows: 

h t {l} = p[A{<r x ®a x ) + B{a y ®a y ) + C{I®I) + D{I®a z + a z ®I)} 
hi{2} = p[E(I ®I) + F(a x ®a v - a y ® a x )\ 
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where a x ,a y , a z are the Pauli sigma matrices and 

1 1 

A — --ie(v 2 )(l + ks(v 2 )), B = --ie(v 2 )(l - ks(v 2 )), k const. 

C = s(v 2 )e(v 2 ), D = ^ic(v 2 )e(v 2 ), 

E = -id(v 2 )e(v 2 ) 2 , F = ^d(v 2 )e(v 2 ). 

where d and c are the elliptic functions dn and cn respectively. 

4.3 Second Order Conserved Currents 

The second order conserved currents can be obtained by the formula (7): 

£{2ei} = t{2e 2 } = £ {a{a x <g> a z <g> a y - a y <g> a z <g> a x ) 

i 

+ (3(a x ®a y ®I -a y ® a x <g> /)} a,P const. 

and 

i{ei + e 2 } = E{7K ® o* ® ° x ) + CK ® <r* ® <r w ) 
+ ^(cr 21 <g> (7* <g> I + <g> 0^ <g> I) 

- + Ofo* ® 1 ® !)} 7, C, 77 const, 

in agreement with [7]- [9]. 
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